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Abstract 

We study the regularization and renormalization of the Yang-Mills 
theory in the framework of the manifestly invariant formalism, which 
consists of a higher covariant derivative with an infinitely many Pauli- 
Villars fields. Unphysical logarithmic divergence, which is the prob- 
lematic point on the Slavnov's method, does not appear in our scheme, 
and the well-known vale of the renormalization group functions are 
derived. The cancellation mechanism of the quadratic divergence is 
also demonstrated by calculating the vacuum polarization tensor of 
the order of A*^ and A~^. These results are the evidence that our 
method is valid for intrinsically divergent theories and is expected to 
be available for the theory which contains the quantity depending on 
the space-time dimensions, like super symmetric gauge theories. 



1 Introduction 

Manifestly gauge invariant regularization is useful for the gauge invariant 
renormalization, but there is not so many invariant regularization for non- 
Abelian gauge theories. The dimensional regularization is the most popular 
method, but is not manifestly invariant method for the theory which contains 
quantities depending on the space-time dimension. The chiral gauge theory 
is a good example of the theory which is not regularized by the dimensional 
regularization because the symmetry depending on 75 is relevant. Frolov 
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and Slavnov proposed an invariant regularization method for this theory 
introducing an infinite number of Pauh-Villars (P V) fields . 

There is a similar problem in the three-dimensional Chern-Simons (CS) 
gauge theory, where the anti-symmetric symbol e^^p is defined by the space- 
time dimension. For this theory, we developed a parity-invariant PV regu- 
larization introducing an infinite number of PV fields p, ^. Using with the 
higher covariant derivative (HCD) method, we can construct a hybrid regu- 
larization of the HCD and the PV without breaking a parity invariance. By 
this regularization, we can classify the universality class of the CS coupling 
shift in the framework of the hybrid regularization ^. 

Another example of the theory depending on the space-time dimension 
is the supersymmetric (SUSY) gauge theory. For this theory, it is believed 
that a manifestly invariant regularization exists but there is few method 
preserving both the gauge symmetry and supersymmetry Since the PV 
regularization and the HCD regularization are defined independently of the 
space-time dimension, the hybrid regularization is expected to be one of the 
candidates for the invariant regularization for SUSY. It is indeed available 
for the three-dimensional SUSY gauge theory and gives the same value of 
the CS coupling shift which is predicted by Kao et al. in the SUSY Yang- 
Mills- Chern-Simons (YMCS) theory 0. 

When our method is extended to the four-dimensional SUSY gauge the- 
ory, it is not obvious whether our regularization works properly because 
the intrinsically divergence appears in four- dimensions though the three- 
dimensional theory is finite. So we have to confirm that our method is 
available not only the finite theory but also the divergent theory to con- 
struct a regularization scheme base on the hybrid regularization that works 
properly whenever the dimensional regularization fails. For this purpose, 
we consider the regularization and renormalization of the four-dimensional 
Yang-Mills (YM) theory in this paper, which is the simplest model of the 
divergent theory. 

The hybrid regularization was originally proposed by Slavnov in 1970s P, 
^, [TD[. When the HCD method is applied to the non-Abelian gauge theory, 
some one-loop diagrams are left unregularized because the HCD term renders 
the gauge propagator less divergent but vertices more divergent. To regular- 
ize these remaining diagrams, other regulator must be introduced. Slavnov 
employed the PV method for the regulator of the remaining diagrams. 

Though the theory was formally regularized by his regularization, it was 
pointed out that his method leads wrong value of the renormalization group 
(RG) (3- and 7-function for the four-dimensional YM theory [rT|. More- 



over, the unitarity breaking was confirmed when his method is applied to 
the quantum chromodynamics |T2|. These problems were occurred by the 
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unphysical non-local logarithmic radiative corrections from the PV determi- 
nant. Namely, the Slavnov's PV field was not the complete regulator for the 
remaining one- loop divergence. 

To overcome this problem, two modifications have been proposed. One is 
to use the dimensional regularization instead of the Slavnov's PV regulariza- 
tion |13[. This scheme does not lead the unphysical logarithmic corrections, 
but is not suit for our aim to apply to the theory which is not regularized by 
the dimensional regularization. Another proposal is to modify the Slavnov's 
PV fields not to lead the unphysical corrections. Introducing 'gauge-fixing 
parameter' for the PV fields, this task is accomplished with preserving the 
gauge symmetry |[T^ . 

On these two proposals, the extra regularization so-called 'pre-regulator' 
is needed to evaluate the divergence of the diagrams. Since it is inserted as a 
partial regulator, the regularization might be inconsistent in the scheme. If 
we use the dimensional regularization as the pre-regulator in the modified hy- 
brid regularization, the scheme is not suit for the regularization of the theory 
depending on the space-time dimension. The extra regulator is sometimes 
unwanted procedure for our aim so we have to develop a procedure without 
it. 

On the other hand, the four-dimensional YM theory leads up to the 
quadratic divergence. Since the quadratic divergence is not ignored in the 
PV type of regularization though it does not appear in the dimensional reg- 
ularization essentially, we also confirm whether it is canceled or not in the 
hybrid regularization scheme. The cancellation, however, is not verified so 
far, though the logarithmic divergence is shown to agree with the physical 



divergence derived from the other regularization scheme [15 



The reason why the cancellation does not confirmed and is not problem- 
atic in references [0, is because they use the dimensional regularization 
as the 'pre-regulator' in addition to the fact that the quadratic divergence 
has no physical meanings. In their scheme, the quadratic divergence does 
not appear in principle though the HCD term leads the non-trivial contri- 
bution of quadratically divergent diagrams to the effective action. So the 
explicit cancellation of the quadratic divergence must be confirmed in our 
regularization scheme for the complete regularization. 

The organization of this paper is follows. In Section |^, we apply our reg- 
ularization scheme to the four- dimensional YM theory writing the explicit 
form of the action. It is confirmed that all the divergences are regularized by 
HCD terms except some diagrams at one-loop level by the superficial degree 
of divergence. A minor modification for the HCD term is given for the com- 
plete cancellation of the quadratic divergence. To treat the unregularized 
one-loop diagrams, the infinitely many PV fields are introduced in a similar 
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way of the CS gauge theory. We check whether the PV fields really cancel the 
divergence by an explicit calculation of one-loop diagrams in the following 
sections. In Section we give a discussion for the contributions independent 
of A. Though each diagram contains up to the quadratic divergence, almost 
the divergence cancels out in total and only the logarithmic divergence re- 
mains. The renormalization is also given to treat the logarithmic divergence 
following the usual procedure. After that, our regularization scheme works 
properly and gives the correct RG j3- and 7-function in this calculation. We 
also show that any explicit pre-regulator is not necessary to see the cancel- 
lation of the quadratic divergence. In Section ^, the one-loop contributions 
of A~^ order are calculated. It is show that the cancellation mechanism of 
the divergence is the same as one of the A'' order. We also show that the 
cancellation mechanism works only when we use the modified HCD term 
in Section ^ Conclusions and discussions are given in the last section and 
some useful formulae are given in Appendices which are employed in explicit 
calculations of one-loop diagrams in Section || and 



2 The Regularization Method 

We consider the hybrid regularization of the YM theory in this section. The 
hybrid regularization consists of the following two steps. First we introduce 
HCD terms. They improve the behavior of propagators at large momentum, 
rendering the theory less divergent at the cost of irrelevant vertices. The 
theory is reduced to superrenormalizable one which has just a finite number 
of divergent loops. As see later, all the diagrams except one-, two-, three- and 
four-point functions at one-loop level are convergent with a suitable choice of 
these terms. Secondly, we deal with unregularized diagrams by a PV type of 
regularization. Since we are considering the gauge invariant regularization, 
the PV regulator must be constructed gauge invariant form and never lead 
any unphysical divergence. 

The generating functional regularized by the hybrid regularization is writ- 
ten by 

Z = f VA^VbVcVc exp[-SA] JJ det A^- JJ det ^'C„ (1) 

where is an action regularized by HCD terms, det and det ^'Ci are 

PV determinants for the gauge and ghost respectively. In the following, we 
consider the regularization in detail and write the explicit form of the HCD 
terms and PV determinants. 
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In four-dimensional Euclidean space-time, the YM theory is given by the 
action 



S — Sym + Sqp, (2) 



where 



1 
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SYM = -r d^xF^F'^^^ (3) 



(4) 



with the field strength F^^ = d^^Al - d^A^ + gf^^^A^^A^ and the covariant 
derivative D^f = S^-^d^ + gf^^^A^^. Here A^, c", c°' and b"' denote the gauge 
field, ghost, anti-ghost and auxiliary field respectively, is the gauge- fixing 
parameter and f"'^'^ is the structure constant of the gauge group SU(A^). We 
abbreviate the color index in the following discussions. 



2.1 Higher covariant derivative method 

The basic idea of HCD method is to regularize the diagrams by improving the 
convergence of propagators with higher derivative terms. When we choose A 
as a cutoff parameter, the most general form of the HCD action is given by 

^HCD = J^I d'xD-F,^D^F^\ (5) 

Though the original propagator derived from (Q) behaves ~ at large 
momentum p, the large momentum behavior of the transverse part of the 
propagator is modified to ~ p-2-2n g^f^gj^ iy^q insertion of (^), but the longi- 
tudinal part is not. So the whole of the propagator still behaves ~ at 
large momentum p. To improve the convergence of the longitudinal part, we 



introduce a higher derivative term H to the gauge-fixing action as follows |T6 
S^^ = j d^x _ hHd^'A^ + cHd^D'^c. (6) 

if is a dimensionless function of d^/A'^ and must contains a higher term than 
Qn I ensure the large momentum behavior of the propagator. 

The most distinctive point of our regularized action is that the higher 
derivative term for the ghost is introduced in @, which is necessary for the 
BRST invariance on our method. There is another choice of the regularized 



action according to the references |]lT], |13| which does not need any higher 
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derivative term for the ghost. It seems simpler than our method and we 
prefer to use their action, but as we see later, the quadratic divergence is not 
completely canceled among the A-dependent terms [0. We will come back 
to this problem in more detail later in Section [4 .21 . 
So the regularized action in (|l]) is given by 

Sa = Sym + 5'hcd + S^p, (7) 

and invariant under the BRST transformations 

(5b^/x = (-D^c), (5b6 = 0, 6bc = -cxc, Sbc = b. (8) 

Here ip x (p means gf^^'^i/j^cj)'^ and then the BRST operator 6b satisfies the 
usual nilpotency 5| = 0. 

The superficial degree of divergence is calculated at 

10 = 4- 2n{L -1)-Ea-(^ + i) E,, (9) 

where L, and Ec are the number of loops, external line of the gauge and 
external line of the ghost, respectively. For all the diagrams higher than 
two-loop (L > 2), n > 2 always gives negative uj. This means that we may 
remove the higher loops by a suitable choice of n. On the other hand for 
one-loop (L = 1), is not always negative by any n. The most economical 
choice is n = 2 and we adopt it. 

The explicit form of H is determined by the behavior of the gauge prop- 
agator which is obtained as follows: 

7(7TA^(^'^- - P,P.) + ^H^^f.P^- (10) 

The first term has the order of the momentum degree of —6, the second term 
must be the same degree or less to ensure the convergence of the diagrams 
except one-, two-, three- and four-point functions at one-loop level; be- 
haves ~ at large p. When we naively remove the higher derivative terms 
taking the limit of A oo, the propagator must recover familiar one; 
converges to unity. So the simplest form of in momentum space is 

It reduces the propagator to very simple form in which the first and the 
second term of ([l0|) has the same denominator. Especially in Feynman gauge 
(^0 = 1) the propagator is reduced to 

A^ 

p2(p4 + A4)^^'^' ^^^^ 
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which renders diagrams simpler so we choose the Feynman gauge in the 
following sections. 



2.2 Regularization of one-loop divergence 

All the diagrams are regularized except one-, two-, three- and four-point 
functions at one-loop level. We need an extra regularization to regularize 



the remaining one-loop diagrams. It was shown in reference [|I0| that the PV 
regularization does not break gauge invariance at one-loop level and therefore 
can be used to complete the regularization of the theory. 
The PV determinant of (|l]) is introduced such as 

oo 

JJdef^Aj- = JJdet"^A+jdet"^A_j. (13) 

j i=i 

The basic idea of our PV regularization scheme is to regularize the theory 
by the pair of two PV determinants det~^~A+j and det~^~A„j. These 
determinants are the generalization of ones in three dimensional CS gauge 
theory: where a pair is needed to make a parity-invariant PV regulator, 
because a part of the action breaks the invariance . The same situation also 
arises in the regularization of the chiral gauge theory |l|, |^ , the regularization 
by the pair is regarded as a more general method than the usual PV regulator. 
So the YM theory ought to be regularized by the similar way of the pair. 

An infinite number is needed to use such PV pairs. Introducing one pair 
corresponds to subtracting double the divergence, so we have to remedy the 
over subtraction by an insertion of another pair of opposite statistics. Then, 
to remedy the over addition we have to introduce the third pair. Such steps 
are repeated alternately until the divergence is removed. The divergence does 
not converged by finite steps but infinite, so we must introduce an infinite 
number of the PV pairs. 

Such steps correspond to introducing fermionic PV fields [aj = — 1) and 
bosonic PV fields {aj = +1) alternately. Following references 0, ^, we take 
the PV conditions for the gauge field such as 

M±, = M|j|, «^. = (-l)lil (14) 

Mj denotes the mass parameter of the PV field A^^, and the determinant is 
represented as 

det A, = j VA^^Vb, exp[-5M, - Sj^^], (15) 
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where 



S 



Mi 



-iH 



X] 



5'Sa 



_6A^{x)6A^{y) 



d^x 







AjAy), (16) 



(17) 



We do not take the summation with the index j in (|T6|) and (|13). bj is an 
auxihary field for and C,j a 'gauge-fixing parameter' for the PV field. 
In usual PV regularization, the 'gauge-fixing parameter' for the PV field is 
ordinary chosen in Landau gauge = 0. In this gauge, however, PV de- 
terminants do not converge formally to a constant and then the anomalous 
divergence appears |l3 . This anomalous divergence contributes to the renor- 



malization and gives the wrong RG (3- and 7-function |TT|, |T4|. To resolve 
such a problem, we have to introduce ^^(7^ 0) in ([T7|). 

H = H{D^ / h?) is a higher covariant derivative term which has the same 
effect as H^d"^ / h?) and satisfies the following two conditions. First, it must 
behave ~ to improve the behavior of the propagator at large momentum 
and converges to unity at the limit A — > cxd to recover the usual PV fields. 
Secondly, it must be invariant under the BRST transformations which come 
from the change of the integration variables 0^ (pj + 6{SB<f)j), where 6 is 
an anti-commuting parameter and (pj denotes a PV field, then (5B0j is given 

by 



^B^j> — X c. 



Ssbj = bj X c. 



On these conditions, H{D^/h?) must be the polynomial of the covariant 
derivative and the simplest form is given by 



(19) 



Consequently the PV determinant (|T5]) is constructed to preserve the BRST 
invariance. 

Similarly, the PV determinant for the ghost is written 



Yl det ^"C, = JJ det ''+^C+i det ^-'C. 



(20) 



1=1 



det^'C, 



VciDci exp 



- / d X CiHD„D^Ci - m, dd 



(21) 
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where we do not take the summation with the index i again, q and q are PV 
fields for the ghost and anti-ghost of mass rrii. We can treat these PV fields 
under the PV conditions m±i = m \i\ and 7±i = (-1)1*1. The HCD term for 
the PV field is inserted as the function H to preserve the BRST invariance 
of 

SsCi = CiX c, JbQ = q X c. (22) 

This HCD term is necessary for the cancellation of the quadratic divergence 
when the mass term is simply introduced in the usual form like m^CjCj, as we 
see in Section ^. 

In the following sections, we confirm that such PV fields completely reg- 
ularize the theory by an explicit evaluation of the one-loop diagrams. In this 
calculation, one of the most important procedure is to summate an infinite 
number of the PV diagrams and derive a convergent function from them. 
Since this procedure is carried out before the momentum integration, all the 
parameters independent of the indices i and j must be chosen to be the 
same. So we assign the same momentum parameter to the internal line of 
each diagram when all the one-loop diagrams are drawn in the next section. 



2.3 Feynman rules 

The regularized action is decomposed into the kinetic part K and the vertex 



part V as follows [T]_ 

d^x^(x)(ir + y + M2)$(x), (23) 



where \E'(a;) and $(a;) denote arbitrary fields and M their mass parameter. 
Since K and V consist of the original part from the YM term (denoting with 
suffix '0') and A-dependent part from the HCD term (with suffix 'A'), they 
are decomposed into 

K = Ko + ^i^A, V = Vo + ^l-A. (24) 

Under this decomposition, the propagators are written in the form 

^ ^ ^1- ,^^^,, A-^ + 0(A-«)^ . (25) 



So the Feynman rules are written by the order of A^^, the quantum cor- 
rections are calculated order by order. The Feynman rules are listed in 
Appendix |^. 
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3 One-Loop Contributions Independent of A 



Now we check whether the PV fields cancel the unregularized divergence by 
an explicit calculation of the one-loop contributions. First we calculate the 
contribution which does not depend on A. 

All the diagrams are manipulated under the following three rules. 

1. Take the same assignment for the internal momentum among graphi- 
cally the same form. 

2. Take the infinite sum of PV diagrams under the PV conditions (|T^) 
adding a 'virtual' PV diagram constructed by taking the massless limit 
of the PV field. 

3. Subtract the 'virtual' diagram to ensure the total contribution as a 
diagram from original fields 

As we mentioned in the previous section, Rule 1 is necessary to find a con- 
vergent function easily after the infinite sum of the diagrams. Rule 2 and 3 
means to divide the contribution of the PV determinants into two parts, the 
part of the infinite sum and of the massless term. 

In the 'virtual' PV diagram, the 'zeroth' field like 'Aq^' runs as the internal 
propagator. We call such a diagram as 'zeroth' diagram in the below. The 
same contribution must be subtracted to maintain the total contribution. 
This procedure is realized by the following equation: 

oo oo 

Yl det ~^A+j det ~^A_j = JJ det A^/det Aq, (26) 

j=l j=-oo 

and so on. The series is defined to give an uniquely convergent function when 
we go to the r.h.s. of this equation. 

In the below, we denote the contribution from the infinite sum by 'mas- 
sive' contribution and the other by 'massless' contribution. The contribution 
from the diagrams of original fields and of massless zeroth PV fields belongs 
to the latter. This decomposition of the contribution clarifies the cancellation 
of the quadratic divergence as we see in the below. 

3.1 Vacuum polarization tensor 

All the diagrams contributing to the vacuum polarization tensor at A° order 
are listed in Figure 0. We denote the quantum correction from the diagram 

(b) 

(b) containing the internal field Aj^ by nj^^(p) and so on. The Feynman 
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(a) (b) (c) (d) 




(h) (i) (j) (k) 



Figure 1: All the diagrams contribute to the vacuum polarization tensor at 
A° order. The wavy line means the gauge field A, the curly line the PV 
field Aj, and the straight line the auxiliary field bj. The rough- and fine- 
dotted line mean ghost c and PV for ghost q, respectively. We use the same 
assignment of internal momenta, where q = k — p. 



gauge ^0 = 1 is chosen for its simplicity in which the propagator is reduced 
to the simplest form as mentioned in the last section. The gauge fixing 
parameter for the PV field is taken as the same value, C,j = 1, for the same 
reason. All the diagrams are divided into three groups and the quantum 
corrections are calculated in each group. 

3.1.1 Gauge field type (A-type) diagrams 

The diagrams (a), (b), (c) and (d) in Figure |l] belong to this group. First we 
consider the diagram (b) that contains only Aj^ field in the internal line. To 
take the infinite sum following (p6D , the diagram of 'Aq/^' must be introduced 
to our calculation. Notice the diagram (a) has the same structure with 

(a) 

the zeroth diagram of (b), both the diagrams are identifiable as Il^,y{p) = 

(b) 

no^i/(p). The infinite sum is taken without any extra virtual diagram. All 
the contribution from these diagrams are written in the form of the infinite 
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sum which is denoted by 'massive' contribution as follows: 



(a) ^ (b) ^ 



9^ 
327r2 



'3M 



.(b) 



j=-oo 



AO 



10 ^^^-+iy^n2MJ-Y-9M^J^'^^^ 



— n — I ^ p^6.^ + O (M-^) 

3 \2MJ 9 360My^ ^ V ; 



49 377r V \ 



(27) 



where we use the PV conditions (0) and the formula (0). Under these 
conditions all the summations are calculated in the same manner with the 
matter field in the chiral gauge theory |1|, ^ using the formula Yl'jL-oo A^+j^ ~ 
. . \ — T. Then the massive contribution converges in the finite M as is shown 

A sm h TV A ° 

in ( P^D in Appendix |C[ The first term of the second line, multiplied by 
6*2, expresses the quadratic divergence at M — > cx). C2 is a dimensionless 
constant which originates from the integration with X = p'^/M'^. 

In the same way, we calculate the corrections from the diagrams (c) and 

(c) (d) 

(d) identifying Ilf,„{p) = Uo^^{p), 



(c) - (d) 



M2 



00 2 
327r^ 5 



+ 



O (M-^) , 



(28) 



This contribution of (ESI) does not lead the logarithmic divergence, but give 
the quadratic divergence as well as (p^). 



3.1.2 Auxiliary field type (B-type) diagrams 

For the diagrams (h), (i), (j) and (k), the zeroth diagram is needed to take 
the infinite sum of these diagrams. There is, however, no diagram which 
is identified as such a zeroth diagram for this group. So we have to add 
the virtual diagrams for (h), (i), (j) and (k) where Bq or Aq field runs as 
an internal propagator instead of bj or Aj. Then the contribution from the 
diagram (h) is calculated 

^ (h) ^ .(h) (h) 

J2 (^j^jM = E i-'^y^jM - noM.(p)- (29) 

j=—oo j=—oo 

The second term of the r.h.s. is a counter term for the zeroth diagram of 
(h) introduced to extract the massive contribution of the first term. Such 
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a counter term is classified into 'massless' contribution. We take the same 
care for tfie otlier diagrams (i), (j) and (k). Q Tlien tlie total of massive 
contributions is 



j=-oo 



. ^ . /(h) (i) (j) (k) \ 

J2 (-1)' ( n,^.(p) + ii,,,{p) + u,,,{p) + n,,,ip)\ 

— — rv-i ^ 



and the massless contributions, 



(h) (i) (j) (k) 



AO 



d^k 1 



(30) 



(27r)4 Pg2 
Here we use (BO) to get the r.h.s. of (I 



{k^ + - 2p^) - k^Qy - k^ky - q^q^ + 2p^py 



(31) 



3.1.3 Ghost field type (C-type) diagrams 

The diagrams containing the ghosts or the PV for the ghosts are classified 
into this group. Since there are some differences in the vertex functions 
between the ghost and the PV field, c and c do not play the role of Cq and Cq. 
The diagram (e) cannot be identified as the zeroth diagram of (f) though the 
both diagrams are the same in graphically. For the diagram (g) , the situation 
is the same as B-type diagrams. So we take the same care in the similar way 
as (|29|) to extract the massive contribution. Then the massive contribution 
of this group is 

(-i)M n,^,(p) + n,^,(p) 

= f 1 In f I^) - 1^ - (p\^ - p^p\ + O im'^) , (32) 

327r2V3 \2mJ 9 90m2y ^ ^ ^ ^ ' 

(k) 

The compensation term for the diagram (k), nop,y(p), is zero because the mass pa- 
rameter Mj is multiphed in the numerator of the integrand all over when we write the 
contribution explicitly in the integral form. But we formally write it in the following 
discussions. 
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and the massless contribution is 



(e) 

n 



ip) 



(f) 



(27r)4 A;2g2 



(33) 



Both the massive contributions from (f ) and (g) give the quadratic divergence 
proportional to C2 after the use of (|86D . These contributions are canceled 
out because they have the same value except the sign, and then the quadratic 
divergence does not appear in (|32|). 



3.1.4 Total contribution 

We can get the total contribution of the vacuum-polarization tensor adding 
all the contributions calculated above. All the quadratic divergence is can- 
celed out exactly and then only the logarithmic divergence remains as follows: 



Total 

n 



AO 



3 \2MJ 9 20M2 J ^'^ ^^^'^^ ^ ' 



+ 



g^ 

87r2 



(pV-PmP-) (Inp' + Ci), (34) 



where the first line comes from the infinite sum (^), (28) and (3C) and 
the second line from the infinite sum (p2D. We use (|7^ to derive the third 
line which is from the massless term (|3ll) and (|33D . As we see in ([76| ) it is 
estimated to have a logarithmic divergence such as Ci = — InA^^ where Ai 
is a parameter of the dimension of one. If we choose M = m = A4, the 
divergent part of the vacuum polarization tensor is calculated 



Total 

n 



iP) 



AO 



div 



g'^c^ 5 



(35) 



This result corresponds to the usual logarithmic divergence of the YM theory 
at the Feynman gauge. 

Here we discuss how to cancel the quadratic divergence in the above cal- 
culation analyzing the contribution from each group of diagrams. We have 
divided the diagrams into the three types, A, B and C. The contribution 
from the massless part does not appear in A-type because the gauge field is 
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identified as the zeroth PV field and all the contribution from this type is cal- 
culated as the massive contribution. This situation is the same as the matter 
field in the chiral gauge theory. On the other hand, b, c and c field is not 
identified as 60, cq and cq respectively, not only the massive contribution but 
the massless one arises in B- and C-type. In C-type, however, the quadratic 
divergence from the massive contribution is cancelled after the summation of 
the infinite diagrams and then the quadratic divergence only comes from the 
massless one. All the arising aspects of the quadratic divergence are listed 
in Table |l] of Section 4. 

The cancellation of the massive contribution occurs between A- and B- 
type. Notice that B-type will be included in A-type after b and bj are inte- 
grated out from the theory, we see this cancellation is essentially the same as 
the cancellation of massive contribution in C-type: the cancellation occurs 
in the same group and the quadratic divergence does not appear in outside. 

On the other hand, the massless contribution is cancelled between B- and 
C-type. We remember that all the massless contribution essentially comes 
from lacking terms to take the infinite sum of the PV diagram, this cancel- 
lation means that the lacking term does not give the quadratic divergence in 
total although each diagram does. 

3.2 Ghost self-energy and vertex correction 

In the above subsection, we show that the vacuum polarization tensor is reg- 
ularized except the logarithmic divergence. Our next task is to treat this 
divergence with the renormalization procedure. The simplest way to renor- 
malize the YM theory is to calculate the contributions of the ghost self-energy 
^(p) and of the gauge-ghost-ghost vertex (TAcc)fi{p)- In this subsection, we 
calculate these contributions at one-loop level to give a renormalization in 
the following. 



We cite the diagram of the ghost self-energy in Figure Since there is 
no contribution from the PV fields we can easily calculate the diagram using 




k 



Figure 2: The ghost self energy loop 
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([7^) as follows: 



AO 



div 



d'^k kp 
(27r)4 A;2g2 



(36) 



In the same way, the one-loop corrections to the gauge-ghost-ghost vertex in 
Figure ^ are calculated as follows: 




/ 



J/\/\/\y\^ 

(b) 



Figure 3: The one loop correction to (TAcc)^l■ 



A" • 3 

(rAc.),(p) =--^p,\nM\ (37) 
div 327r2 

Here we use (0) again. 

3.3 Renormalization 

Now we give a renormalization procedure to absorb the logarithmic diver- 
gence which we calculate in the preceding subsections. Since we are only 
considering the one-loop corrections, the usual renormalization procedure of 
the YM theory can be used here. In that procedure, it is well known that all 
the divergences are renormalized by the three renormalization constants Zi, 
Z3 and Zc as follows: 

11 i_ 

_1 _3 l'^^) 

^bare = ^3 ^ ^) Qhare = ^1^3 ^9, ^Obare = ^3^0- 

Here we denote the bare parameters with the index 'bare'. Expanding the 
renormalization constants with h such as Zi = ff'zf''^ = 1 -|- hzf''' + 0(/i^), 
the first order of h corresponds to the one-loop corrections. Then we get the 
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following equations: 



A" 

div 
A" 

div 
AO 

div 



(39a) 
(39b) 
(39c) 



Comparing these equations with the results (|35|), (p^) and (^), the renor- 
malization constants are easily calculated as follows: 



1 l«o=i I67r2 3 ' ^ l«o=i 



167r2 3 ' " 



Co=i I67r2 



In M. 

(40) 



From these results, the j3- and 7-function at a renormalization point are 
calculated 

11 

i&=i i67r2 3 



/5(^7,eo) 
7a (^,^0) 



/i (9 In 
2 (9/i 



■c. + 0(/), 



5 



167723 



77C. + 0(/). 



(41a) 
(41b) 



These values accord with the familiar value at Feynman gauge |18, 19, 20 1. 



4 Contributions from HCD Terms 

In the last section, we have seen that the quadratic divergence is completely 
canceled and the correct RG (3- and 7-function are given with our regular- 
ization scheme by the calculation of the one-loop contributions in A" order. 
In this section, we mainly consider the contribution from A-dependent part. 

We remind that the HCD terms render the renormalizable theory into a 
super-renormalizable one, it is reasonable to consider that all the divergence 
from higher than two-loop level in the renormalizable theory is translated 
into some divergence at one-loop level in the super-renormalizable one. So 
there is some divergence in A-dependent terms and it must be regularized by 
our regulators. We now check such divergence is regularized at least in A~^ 
order. 



4.1 Vacuum polarization tensor in A ^ order 

A new diagram listed in Figure ^ arises in this order in addition to the 
diagrams in Figure 0. This is from the irrelevant vertex which maintains the 
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(1) 

Figure 4: Diagram appears in the order. 



BRST invariance. We classify this diagram into B-type because it contains bj 
field in the internal line. The contribution is calculated in a similar way with 
the case of A° order under the rules mentioned on the top of Section ^. Since 
the contribution is divided into two parts following the calculation rules, one 
is the infinite sums and the other is the massless terms, all our tasks are 
concentrated on the calculations of the infinite sums (arising from the PV 
fields like (|30D) and of the massless terms (from the counter terms such as 

For the contribution of the infinite sums, all the terms are written by the 
following formula: 



E 



A4 J (27r)4 {k^ + Mlf{q^ + Mp^ 



i_(^M2[3-f]+o(M^[^-i])), (42) 



where r and s denote the order of and in the integrand, and run in the 
region < r < 5 and < s < 10 — 2r respectively. We use the computation 
rules discussed in Appendix to get the second line and [ ] denotes the 
Gauss' notation. Its contribution depends on the order of the limit M and 
A, but the difference does not affect renormalization group functions |Tl|, we 
may restrict our calculations within the ratio M/A = constant. Then (|42|) 
behaves 

~ M^t^^t] +0 (^M^t"^]) , (43) 

and only the terms having s = give the divergence in the limit M — oo, 
they are calculated in (P7|). 

Similarly the contributions from the massless terms are written in the 
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following form using the formula ([75| ) and (|7^), 



1 r d'^k ■■■k 



Ml 
1 



M 



2 3- 



+ P 



2[3- 



(44) 



where < s < 6. Also in this case, only s = terms give the divergence 
under the condition M/A = constant if we identify M. with M. 

Both in (|42D and (0), the logarithmic divergence appears in the form 
of InM^/A^. Since it vanishes in the ratio M/A = constant there is no 
logarithmic divergence and only the terms containing no external momentum 
p^ in the numerator give the quadratic divergence in this order. So we only 
consider s = terms and confirm whether the divergence cancels or not by 
calculating each diagrams. 

Now we calculate the divergent contributions in the below. For the di- 
agrams in A-type, since the diagram (a) and (c) play the massless diagram 
respectively we take the infinite sum without any extra diagram as in the A*^ 
order. The divergent contributions from (a), (b), (c) and (d) are 



E 



(d) 



9 



div 



87r2 154 



(45) 



Here we use (|^) and C4 is a dimensionless constant. 

For the diagrams in C- and B-type, we take the infinite sum adding the 
external diagrams as in 



(<=) ^ . /(f) (g) \ 

i=—oo ^ 



div 



d^k 



A4 J (27r)4A;2g^ 



{2k^6^, - 4k%k,) , (46) 



j=—oo 



(h) (i) (j) (k) (1) 

nj>v(p) + ^j,,u{p) + ^j,iu{p) + nj^^(p) + nj^j,(p) 



div 



^^c^ 9 
87r2 154 



A4 



d^A; 



(27r)4 A;2g2 



{2e5^,-Ak%k,). (47) 



The first term of the r.h.s. in 



arises from the infinite sums of PV 



diagrams, and the second term from the massless terms. In 



we add 
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Diagram type 


A° order 


A ^ order 


massive massless 


massive massless 


„ A- type 
Gauge 

B-type 


+^M^C2 


-^M6C4 


Ghost C-type 


-IM^ 


-ToIa^-^' 



Table 1: Quadratic divergence at A° and A ^ order, 'massive' and 'massless' 
means the contribution from the infinite sum such as (|30| ) and the massless 

2 

terms such as ( |3TD respectively. The coefficient ^^5^^ are abbreviated from 
explicit values. 



the contribution from the diagram (e). The each diagram (f) and (g) gives 
the quadratic divergence proportional to the constant C4 with the usage of 
(|87|) , but they have the opposite sign and cancel each other. Then only the 
quadratic divergence from the massless terms remains. This situation is the 
same as we see in (p3D and also at the A*^ order. 

All the quadratic divergence of the order of hP and A~^ are listed in 
Table |l|. It is easy to see from p5|), (|46| ) and p7|) that all the divergence of 
the vacuum polarization tensor at A~^ order disappears from the theory. It 
clearly shows that the divergence at this order is removed completely. This 
cancellation mechanism is exactly the same as we observed in Section ^ B- 
type contribution gives the counter terms for A- and C-type. Considering 
A-type diagrams do not give any massless contribution in all orders of A, the 
same mechanism is strongly expected to work in the higher order of A~^, e.g. 
in A~® order. 

Such a cancellation mechanism is observed without an explicit identifica- 
tion of the parameters M = m = M. which is necessary for the renormaliza- 
tion, because the cancellation is occurred between the divergence described 
by the same parameters: massive divergence is canceled by massive one and 
massless divergence by massless one. 



4.2 Inconsistent higher derivative action 

There is another choice of gauge-fixing action but @. Remembering that the 
higher derivative function for the gauge-fixing action is added to improve the 
convergence of the longitudinal part of the gauge propagator, we can choose 



the action with a higher derivative function / according to the references [11 
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|T3| 



as follows; 



4p = j b^b - + cd^D'^c. (48) 

This action is extended to the Pauli-Villars field using a higher covariant 
derivative function /, 



2 



(49) 



These actions give the same propagators for the gauge field and PV for gauge 
as we get with the actions (^ and (p!7D under the condition f = H and f = H, 
but change some propagators and vertices. Following points are modified by 
the usage of the action ( ^8]) and (^); 

1. The A-dependent parts of propagators (Ajbj) and (bjbj) and vertices 
(AbjAj) and (AAbjAj) are modified. Especially, these vertices do not 
give any A-dependent term; 

2. New vertices (Abjbj) and (AAbjbj) appear because of /. These vertices 
give the new diagrams. 

3. The ghost and its PV field have no A-dependence. 

From the first fact, (^Tj), which is the total contribution of the diagrams 
(h), (i), (j), (k) and (1), is changed. The diagram (1), especially, does not 
appear because of the absence of the vertex {AAbjAj). Recalculating these 
diagrams with modified propagators and vertices, we get 

^ . /(h) (i) (j) (k) 



j=—oo 



div 



a'^Cv 6 , . 



87r2 154 



The second fact leads the new diagrams listed in Figure ^ Since the 
diagram (m) only gives the higher order of A~^ and the diagram (n) does 
not give any quadratic divergence, the diagrams we have to calculate here 
are the three diagrams of (o), (p) and (q). Each diagram gives the quadratic 
divergence after the infinite sum but all of them are canceled in total, the 
massive contribution does not appear from these new diagrams. Not only 
the massive contribution but massless one is not arising from these diagrams: 
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(m) (n) 




Figure 5: Diagrams generated by the new vertices appearing in (^8[) and 
PS|). These diagrams only contribute to the A depending terms but do not 
give any quadratic divergence to in totaL 



all these diagrams contain mass parameter Mj in the numerator of integrand 
so their zeroth diagrams vanish taking the massless limit of them, like the 
diagram (k) in Section |^. 

Since the ghost and its PV do not give any A-dependent contribution in 
this case, the quadratic divergence in the A~^ order only comes from (^51) and 
(|50|) . Comparing these two contributions, we find the quadratic divergence 
does not cancel in this order. 

The reason why the action (^) and (^9]) fail the cancellation of the 
quadratic divergence in A~^ comes from the fact that the function / gives 
a different effect to the mass term of the PV field for ghost, with the case 
of H. Certainly, ((|9])) is translated into (D (([l^)) by the redefini- 
tion of c — > cH (cj — ^ CiH) and h hH {bj —>■ bjH) under the condition 
f = H {f = H), but the PV for ghost does not: the mass term gets the 
A-dependence and def^'Cj does not coincide with (pT]). This fact says that 
as long as we construct PV fields by simply adding an usual mass term, (HSD 
and (^9]) are not the correct higher derivative regulator. If we want to use 
these regulators, we must add a A-dependence to the mass term of the ghost 
PV field. 

From an another point of view, it is simply recognized that our actions 
are natural. The main purpose to introduce a higher (covariant) derivative 
function to the gauge-fixing action is to improve the convergence of the prop- 
agator of the gauge field as we explain in Section P?T] . So the function must 
regularize the gauge field. The function / (/) surely improves the conver- 
gence of the propagator by regularizing the gauge-fixing parameter (^j); 
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but it is not the regulator for the gauge (PV) field. On the other hand, the 
higher (covariant) derivative function H {H) regularizes i^jfi) by multi- 
plying to the cross term with the auxiliary field b {bj) and gives the complete 
cancellation of the quadratic divergence without any A-dependence in the 
mass term of the ghost PV field. 



5 Conclusion and Discussion 

In this paper, we check the consistency of the hybrid regularization in the 
four-dimensional Yang-Mills theory when we use the regularization scheme 
consists of the higher covariant derivative term and an infinitely many Pauli- 
Villars fields that we applied to the three-dimensional Chern-Simons gauge 
theory. By an explicit calculation of the vacuum polarization tensor, we get 
the correct factors of the renormalization group f3- and 7-function. Further- 
more, the cancellation of the quadratic divergence is also demonstrated in our 
method with an expansion of A~^. These facts show that our regularization 
method is available for a divergent theory not only for a finite theory. 

In our calculation, all the diagrams are classified into the three types and 
their quantum corrections are calculated separately in each type. In that 
calculation, all the contributions are separated into 'massive' and 'massless' 
contributions along the manipulation rules given in Section 0. As a result, we 
can clearly confirm the cancellation of the quadratic divergence: the 'massive' 
contribution is cancelled between A- and B-type diagrams and the 'massless' 
between B- and C-type. This clear cancellation mechanism works identically 
in both the order of A° and A~^ as we listed in Table |I]. 

Such a mechanism is expected to occur in the quadratic divergence in the 
higher order of A~^, e.g. A~^, because the massless contribution of A-type is 
always absent. This argument is proved if it is shown that C-type diagrams 
do not give the massive contribution in any order of A^^. 

The higher (covariant) derivative functions for gauge-fixing terms, H and 
H, especially, are important for the cancellation of A-depending contribu- 
tions. The actions ( ^HD and ( ^91) sure give the cancellation in A*^ order, but 



do not in A~^. The reason is that the mass term of the ghost PV does not 
accord with ours under the redefinition of fields. In other words, the func- 
tions / and / are only the higher (covariant) derivative for the gauge-fixing 
parameters and they do not regularize the gauge and its PV fields completely. 
Though the higher (covariant) derivative terms are originally introduced to 
improve the convergence of the gauge and PV field, such functions conse- 
quently must be inserted to regularize the gauge and PV field like the action 
(|^) and (|T^, as far as we introduce an usual mass term like mfcjCj. 
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In usual case, some extra regularization so-called 'pre-regulator' is intro- 
duced to compare all the divergence, in our method, however, such a regu- 
larization is not needed. This is because that our method contains a scheme 
corresponding to a pre-regulator. In our calculation, all the divergence ap- 
pears in the constants of integration which is derived from the differential 
equations of X like (|69|) or (^). As a result, 'massive' divergence is char- 
acterized by M or m and 'massless' by A4, and then we get the correct 
renormalization supposing M = m = A4. This assumption will corresponds 
to the conventional pre-regulator in reference where it is only needed to 
give a rigorous arguments. So our method is an alternative procedure which 
does not need any pre-regularization. 

The very reason that we want to avoid the pre-regulator is that the proce- 
dure may break the invariance which we would like to preserve in the theory. 
But in our method, since a scheme corresponding to a pre-regularization is 
in the usual regularization method, we do not worry about the symmetry 
breaking. 

Finally we comment on the problem of the overlapping divergences; how 
to treat the one-loop divergence with external P V fields ||2l|, ^ . It is difficult 
to remove the divergence with our present method, however, the PV pair 
such as (|13D will give the key to the problem. Since both the fields give 
the same diagram with each external field, if we can find a PV pair whose 
diagrams cancel each other like the parity-odd contributions in reference 
where the pair is connected by the parity-transformation, all the one-loop 
divergence with external PV fields are canceled and the problem will be 
overcome. Unfortunately, such a PV pair is not found up to the present. 
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A Feynman Rules 
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A.l Propagators 

A. 1.1 Propagators for original fields 



Xab 



A' 



y(p4 + A4) 



9^,u + (^0 - 1) 



p2 



1 - 



p 
A4 



(51) 



A = !^ ^ !^ (i 



p^ 



p 

2M 



(52) 



-6^" -6' 



ah 



p2jj p2 



1 - 



p 

2M 



(53) 



A. 1.2 Propagators for PV fields 



A^[p-igt^^ -pt^p'^) 



+ 0- 



p2 



2^ / 1 ^6 



p2 + M2 



+ 



A4 p2 + M| 

1 



p2 + e,M| 



A^p^ + ^jMf 



(54) 



4 .'^ 



p2i72 + ^ . JVf 2 p2 + ^^.JVf 2 1^ 2A4 p2 + ^ ) 



p' p^ 



(55) 



-5^^M] 1 



pm^ + e^Ml p2 + ^^.JVf 2 \^ A4 p2 + ^ .^^^2 



P" 



(56) 



„a P -b 
Cj Cj 



-5° 



-(5 



p2i7 + mf p2 _|_ 77^2 Y 2A^ p2 _|_ ^2 



1 - 



P" 



2 / ' 



(57) 
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A. 2 Vertices 

A. 2.1 Three-point vertices 




AO-S A a2 




A "3 



^9 fa\a2ai 



A4 



}9_ jaia2a-j, 



A4 



sym 



- {pIpi - ptP3)^.29^3^.l - {pIps - P2P2)mi^m2M3 - (P2P2 - ptpi),.^9f.^ 



^2 



+ {pl+pl)iP3-Pl)^.2iP^^^lPl^^3 - PsPigf^i^^i) 

+ (Ps + P2)(P2 - Pa) Ml (P2M3P3A.2 - P2P3^/.3/^2) 

+ (P2 +Pl)(Pl -P2)m3(P1M2P2^.i -PlP2^M2w) ' (58) 



where 



J sym 



imphes that symmetrize aU the indices. 




0213 



^MiM3 + ^ (P2^MiA.3 + (pI + PDiPS - ^2)^1^3^.3) 



(59) 



411 



P3 



(60) 



^^2 



r.a3 
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Ml 



P2. 



1 + 



2A4 



(Ps + P3P2 + P2) 



(61) 



A. 2. 2 Four-point vertices 




4^4 /lai 




) + {Pl+P2) (c/miM3£/m2/^4-^M1M4£/m2w) 



+ 8(^1+^2)^1 {P'^^,■AP2Pi9^l2^^^ - P2^hPa^i2) - P3m4 (P2P35'm2/^3 

+ {fl'M2M3(P4MiPlAi4 + PlP4fl'/.iM4) - fi'M2M4(P3MiPlAi3 + PlP3fi'MiM3)} 

- 4p2(2pi +p2)^2(^MiM4PlM3 - ^MiM3PlM4) 
+ 2pi^i {P3fis{P2P'igf,2fi4 -P2M4P4M2) -P4M4(P2P3fi'M2M3 -P2/.3P3M2)} 

+ 4(pi + pa)^ {^m2M4P4mi (Ps + 2P4)m3 - ^m2M3P3mi (P4 + 2^3)^4 } , (62) 

J sym 




^AX4 



P2 P3 



= g'r'^-'r^^' g,,,, + ^ { ^2 + pIpI + ^3)^7^.4 

+ {pI + {P1+P2Y +Pl)ij>2+P3-Pi)t,^ij>2+P3+Pl)t,}^ 



+ (1^4) (63) 
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4 "4 



Pi 



"^2 



A, 



0-3 



^9 ^a\a2b ja^a^b 



2A4 



114 



+ {P3 +Pa- P2)f,i (P3 + Pi - P2)f,4P3fis 



+ (1^4). (64) 



Where (1^4) means the same expression exchanged all the index 1 with 4. 



B Momentum Integration 

All the momentum integrals that we consider in this paper are reduced to 
the general form of 



k/ii ■ ■ ■ 



(65) 



where q = k —p and p is the external momentum. This expression is written 
as follows using the Feynman parameterization and introducing z which is a 
source of the momentum: 



d^k 6 6 



(27r)4 6z^, 6z^^ Jo 



N 



i=l 



exp 



N 



N N 



N 



-k^ ^(a, + A) + k{^P ^ A + z) - Y,i(^i + -P^J2f^' 



1=1 



1=1 



i=l 



1=1 



z=0 

(66) 



First we take the momentum integration with k using the gaussian and the 
differential with z. Inserting an unity 1 = dp6 {p — J2 '^i ~ J2 f^i) ^ 
rescaling the parameters such as aj — > pa^ and Pi pPi, the expression is 
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written as follows: for odd n — 2m — 1 

1 



'■Ml'"M2n 



E 

k=l 



24+TO— fc^2 



sym 



2fc-l 



X dpp' exp J-p ^ Qij ^ /3j + ^(oii + Pi)x,i^ 
and for even n = 2m 

1 



(67a) 



fc=0 



24+m-A;^2 



' PlJ-2k^lJ-2k+llJ-2k+2 ' ' ' '^M2m-l/i2m] 



sym 



X 



dpp^ exp [-p ^ ttj ^ /?j + ^(ttj + I 



(67b) 



J sym 



as 



where i runs from to N. We define / = 2N + k — m — 3 and [ 
the symmetrization about indices hke [Si_,i^^Pu,^]^y^ = 5^^^^p^^ + 5^^^,^p^,^ + 
^ixzij.iPix2- We take ^ ctj + ^ /3j = 1 because of the (5-function. 

To calculate the expression, we consider the p-intcgration previous to the 
parameter integration of and f3i. First wc consider the case when / is 
positive value. With an integration by parts we get, 



poo 

I dpp^ exp [-pb{aX + Y}] 
Jo 



p' 



b{aX + Y} 



^-pb{aX+Y} 



OO /.c 

+ / 

Jo 



dp 



IP 



h{aX + Y) 



^-pb{aX+Y} 



(68) 



The first term of r.h.s. vanishes under the condition h{aX + y} > 0. Inte- 
grating by parts recursively, we get the following result: 



-1 



.1 Ql 



dpp'exp[-p6{aX + r}] = i^ 



{dxy 



aX + Y 



(69) 



Going to the last expression, we use the formula 



/■OO 

/ dp exp [-ph{aX + Y}] 
Jo 



1 



b{aX + Y}' 



(70) 
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For negative / using the relation 

and (0), we get a differential equation about X. The differential equation is 
solved by integrating with X recursively, and the solution is formally written 
by 



Here Cj arises from the z-th integral and represented using the parameter X 
which has the same dimension with X as 

Ci = -lnA', C2 = X, ^3 = -^, ^4 = ^, •••• (73) 

Using the results, we calculate the easiest formulae of = 1 and Xi = 
taking X = as follows: 



m = "1^(1^ + Ci), /.(O) = - J^(l^^^' + Ci). (74) 



We can calculate the more general n of this case 

^Ml-Mn (0) ~ (^[l+f ] + P (^[f ] + ■ ■ ■ + Jnji^-^ + CO'^St. 



(75) 



where [ ] denotes the Gauss' notation and Cj is represented using the pa- 
rameter A4 of mass dimension of one by 

C, = -lnM\ C2 = M\ C3 = — -, C4 = --, (76) 

4 io 

The parameter M. goes to oo just like the cut-off parameter. Under this 
parametrization, ( ffil) is read as the dimensionless divergence of ln(A^/p). 



C Infinite Sum 

Here we give a formula that we use the calculation of an infinite sum of PV 
fields. Under the Feynman gauge, all the expressions of the diagrams listed 
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in Figures |1| and | are reduced to the following forms; 



oo „ 



k ■ ■ ■ k 



j=-oo 



(27r)4(A;2 + M|)^(g2 + M|)^' 



(77) 



At the beginning, we consider r = case. Inserting (^) to this expression, 
we find that the most important calculation is the summation with j which 
is given after the p-integration. 



E 



j=-oo 



dpp'exp [-pM'^b{aX + |jf }] 



(-1)' 



7r2 7it' 



aX 6 360 



ai(M%y+^ {dxy 

where we take Mj = M\j\ and define a, b and X as follows: 



aX + ■ 



(78) 



i=l 



{-ly 



(79) 



and 



Going to the r.h.s. of (^ we use the formula Y.j& ~ Asinh-nA 

expand it. The negative / means multiple integral of order / as mentioned 
Appendix |B[ In the case of / = —1, for instance, we get 



^ /"OO 1 

^ / ^exp [-pM%{aX + |jf }] 



j=-oo 



_lnX + -aX-— -2^2 
6 720 



a^X^ - aCi 



^0) 



where Ci is a dimension less constant and determined by the relation 



^ tanh |-\/axj such as 



/ aX sinh vrV aX 



Ci = - In — a . 



7r 



For r = 1 the key object is the following equation: 

E / dpp'exp [-pM^H^^ + Ijf }] 



(821 



]=-oo 
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Notice the relation 



d 



dpp'-^exp [-pM'^b{aX + |jf }] 

POO 

= -b{aX + |jf } / dpp' exp [-pM%{aX + |jf }] , 
Jo 



and using ( |78D we get 

OO „QQ 

J2 / dpp' exp [-pM'b{aX + m] 



J = -CO 



9 



9X 



d 



dX 



i-i 



1 77r2 ^ 

\ aX + 

aX 6 360 



Using (^) and (|8^ iteratively, we get the formula for r > 1 

OO 

^ (-lyMf / dpp' exp [-pM'b{aX + |jf }] 



j=-oo 



\l+r 



q=0 



:X 



r—q 



qj il-q)\ 

oX ~ y ^ 360 



dX 



aX + 



^5) 



Here we give explicit calculations which we use in Section ^ and For 
A° order, all the integrals are given by (^) with N = 1, we have to calculate 
up to = 2 with r = and n = with r = 2. Since these integrals give the 
p-integration with I = —1 and I = —2 there appear dimension less constants 
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Ci and C2 as follows: 



j=-oo 



^2ln(^)+2 + ^+0{M-^) 



\2MJ 36M2 



, (86a) 



2\ _ PfJ-l 



167r2 



j=-oo 



2^2 



(86b) 



:;=-oo 



+ 



327r2 



M2 



■a 



30p2 3 \2M) 9 360M2 



, (86c) 



00 _ 
E (-1)W|/(M|) = - 



3=-oo 



167r2 



30 6 360M2 ^ > 



(86d) 



Similarly for A * order all the integrals are given by = 2. The most 
divergent parts are calculated 



E (-l)%,..^,o(M?,M|) 



3=-<x> 



j=-oo 



j=-oo 



j=-oo 



2%2 



2772 



C4 + O (M^) , 



E {-iyMp,,...,,{MlMf) 



[^filH2 ' ' ' ^/iy/islsym 3M^ 

2%2 2772 



(87a) 
(87b) 



Y.{-iyMp,,...,,{MlMf)^ 



[^/iHU2 ■ ■ ■ ^/^SMelsym 6M® 

2^^r2 2772 

^3/^4 J svm 

6M6 



C4 + O (M^) 



267r2 2772 



C4 + O (M^) , 



(87c) 

(87d) 
(87e) 



3=-<x> 
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where C4 is the integration constant arising from the fourth order integration 
of X which is caused by the p-integration with / = —4. 
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